COMPLEX ZEROS OF REAL ERGODIC EIGENFUNCTIONS 



STEVE ZELDITCH 

Abstract. We determine the limit distribution (as A oo) of complex zeros for holomor- 
phic continuations (p^ to Grauert tubes of real eigenfunctions of the Laplacian on a real 
analytic compact Riemannian manifold {M,g) with ergodic geodesic flow. If {fj^} is an 
ergodic sequence of eigenfunctions, we prove the weak limit formula j-[Z^c ] —f ;^i9i9|^|g, 

where [Z c ] is the current of integration over the complex zeros and where d is with respect 
to the adapted complex structure of Lempert-Szoke and Guillemin-Stenzel. 



1. Introduction 

A well-known problem in the geometry of Laplace eigenfunctions is to determine the 
asymptotics of the volume of their nodal hypersurfaces in the limit of large eigenvalues. At 
the present time, the best result for compact real analytic Riemannian manifolds (M, g) of 
dimension m is the estimate 

(1) ciX<n'^-\z^j<C2X. 



due to Donnelly- Fefferman jPFj . Further background and references can be found in |DF| l,TH 
INPSoj . In this article we are concerned with the yet more difficult problem of the asymptotic 
distribution of nodal hypersurfaces, i.e. with integrals of continuous functions over nodal 
hypersurfaces (cf. fllO|) ). This problem is much too difficult for real hypersurfaces, but it 
turns out to simplify quite a bit if we complexify the problem, i.e. analytically continue the 
eigenfunctions into the complexification of M. Our main results determine the asymptotic 
distribution of complex nodal hypersurfaces for eigenfunctions of real analytic Riemannian 
manifolds with ergodic geodesic flow, or more generally for any sequence of quantum ergodic 
eigenfunctions. 

To state our results, we need some notation. Let (M, g) be a real analytic Riemannian 
manifold of dimension m, and consider an orthonormal basis of real eigenfunctions 



^gV^j = ^j^j^ {^j,^k)L^(M,dvoig) = / (pj{x)(pk{x)dvolg{x) = 6jk 



M 



of its Laplacian A^. We use the sign convention for which the Laplacian is positive; we 
often write as A when the metric is understood. As reviewed in ^ a real analytic 
manifold possesses a Bruhat- Whitney complexification Mc, that is, a complex manifold in 
which M embeds as a totally real submanifold. This complex manifold may be identified (by 
means of the complexified exponential map) with a ball bundle B*^M inside the cotangent 
bundle T*M, equipped with a complex structure Jg adapted to the metric in the sense of 
Guillemin-Stenzel pTTlirT^ and Lempert-Sz5ke [ESIIE^. We denote by l^l^ = J^ij 
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the norm-squared of ^ G T*M with respect to the metric, and by d the Cauchy-Riemann 
operator with respect to Jg. The maximal ball bundle B*^M to which this complex structure 
extends is known as the Grauert tube of {M,g). The natural symplectic form uj and the 
complex structure Jg endow B*^M with the Kahler metric ujg = j9c)|^|^. 

By a theorem due to Boutet de Monvel |Bouj (see also |(tS2| K^L Sj) the eigenfunctions 
possess analytic continuations y?^ to the maximal Grauert tube. The complex nodal hyper- 
surface of an eigenfunction is defined by 

(2) Z^c = {Ce BIM : ^l{0 = 0}. 

There exists a natural current of integration over the nodal hypersurface in any ball bundle 
B*M with e < eo , given by 

(3) {[Zr],^) = ^[ dd\og\^^,\'A^= [ V, ^GP(--i'™-i)(i?:M). 

^ 27r Js^ j^ Jz c 

^\ 

In the second equality we used the Poincare-Lelong formula. The notation T)^"^~^'"^~^\B* M) 
stands for smooth test [m — 1,171 — l)-forms with support in B*M. 

The nodal hypersurface Z^c also carries a natural volume form \Z^c\ as a complex hyper- 

^m— 1 

surface in a Kahler manifold. By Wirtinger's formula, it equals the restriction of to 
Z^c. Hence, one can regard Z^c as defining the measure 

(4) (|Z^cU) = y^ "^J^^frjy,^ veC{BlM). 

We prefer to state results in terms of the current [Z^pc\ since it carries more information. 
We will say that a sequence {^j^} of L^-normalized eigenfunctions is quantum ergodic if 

(5) (^^^,,^J^_J_|^^^^,,<i^. ^AenM). 

Here, \Ef*(M) denotes the space of pseudo differential operators of order s, and dfj, denotes 
Liouville measure on the unit cosphere bundle S*M of (M, g). More generally, we denote by 
dfir the (surface) Liouville measure on dB*M, defined by 

(6) d/i. = ^ on dBlM. 
Our main result is: 

Theorem 1.1. Let {M,g) be real analytic, and let {fjf.} denote a quantum ergodic sequence 
of eigenfunctions of its Laplacian A. Let {B*^M, J) be the maximal Grauert tube around M 
with complex structure Jg adapted to g. Let e < eo- Then: 

In other words, for any continuous test form ijj G T''('"~^''"~^)(i?*M), we have 



■^i -I Z c ^ J B*M 

3k 
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The limit (1, 1) form ^dd\^\g first arose in [LSHIGbT] for a different reason which is reviewed 
in H2.ll As a corollary we obtain a similar result on the integrals of scalar functions against 
the measures \Z c 1: for any (p G C(B*M), 

1 r , ,m~l ■ p , ,m—l 

\jzc ("^-1)! T^JBtM (m-1)! 

As is well-known, ergodic sequences of density one in the spectrum arise when the geodesic 
flow is ergodic, and an entire orthonormal basis is ergodic when the Laplacian is quantum 
uniquely ergodic |CV| IShn] U\ IZ2j. Thus, we obtain the titled result: 

Corollary 1.2. Let {M,g) be a real analytic with ergodic geodesic flow. Let {(fj^,} denote 
a full density ergodic sequence. Then for all e < cq, 

-^[Z^c ] -> -dd\^\g, weakly m V'^^'^\B*M). 

By the unique quantum ergodicity result of E. Lindenstrauss |Linj . the zero currents 
of the full sequence of Hecke eigenfunctions on arithmetic hyperbolic surfaces satisfy the 
limit formula in Corollary 11.21 The Rudnick-Sarnak conjecture |R,Sj that negatively curved 
compact manifolds are quantum uniquely ergodic would imply that the zero currents of the 
full sequence of eigenfunctions should satisfy the limit formula on such spaces. 

Theorem 11.11 also has implications for a general Riemannian manifold {M,g). The or- 
thonormal basis of a general Riemannian manifold is not quantum ergodic and moreover 
the complex zeros do not generally tend to the limit ^dd\^\g (a fiat torus provides a simple 
example where zeros concentrate on complex hypersurface). However, in a precise sense, a 
random orthonormal basis of L^{M) (adapted to A^) has the quantum ergodic property, and 
hence the complex zeros of the basis functions will satisfy the limit formula of Theorem 11.11 

To state the result, we recall the definition and results on these random orthonormal bases 
from [Z3]. We partition the spectrum of y^A^ into the intervals Ik = [k,k + 1] and denote 
by Ilfc = E{k + 1) — E{k) the spectral projections for y^A^ corresponding to the interval 
Ik- We denote by N{k) the number of eigenvalues in and put Ti.^ = ranllfc (the range 
of Ilfc). Hk consists of linear combinations X]j A e/^ '^j^j eigenfunctions of a/A^ with 

eigenvalues in Ik- We define a random orthonormal basis {Uk'Pj} of Tik by changing the basis 
of A/A-eigenfunctions {(pj} of A in Hk by a random element Uk of the unitary group U{Hk) 
(equipped with its normalized Haar measure du^ of the finite dimensional Hilbert space TY^. 

We then define a random orthonormal basis of L'^{M) (adapted to Ag) by taking the prod- 
uct over all the spectral intervals in our partition. That is, we define the infinite dimensional 
unitary group 

t/(oo) = U^^.UiHk) 

of sequences (f/i, U2, ■ ■ ■ ), with Uk G U(Ti.k), and equip f/(oo) with the product measure 

duoo = ^T=id^k- 

A random orthonormal basis = {{Uk'Pj)} of L'^{M) is thus an orthonormal basis obtained 
by applying a random element U G U{oo) to the orthonormal basis $ = {^Pj} of eigenfunc- 
tions of 
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In |Z3j , it is proved that random orthonormal bases satisfy the following variance asymp- 
totics: 

To be precise, in |Z3j it is assumed that the widths of the intervals Ik increase to infinity, but 
the more recent strong Szego asymptotics of ILR,Sj allow one to prove the same result 
for intervals of bounded width such as the Ik above. By the strong law of large numbers (see 
|Z3j ) it follows that with probability one, a random orthonormal basis of L^(M) is quantum 
ergodic. We thus have: 

Corollary 1.3. Let {M,g) be any real analytic compact Riemannian manifold. Then with 
probability one, a random orthonormal basis {ipj = U(fj} of L'^{M) as defined above satisfies 

-^[Z^c ] -dd\^\g, weakly m V''^^'^\B*M), 

for a full density subsequence {"ipjk}- 

This gives a kind of almost sure improvement of the complexification of Theorem 14.3 of 
Jerison-Lebeau jJL from an inequality to an asymptotic formula. 

1.1. Discussion and outline of the proof. In summary, the complex zeros of a quantum 
ergodic sequence become equidistributed with respect to the (1,1) form ^dd\^\g. As men- 
tioned above, the Kahler form on B*^M associated to g is the (l,l)-form Ug = We 

observe that dd\^\g is singular relative to the Kahler form along the zero section, i.e. the to- 
tally real submanifold M (the geometry will be reviewed in This singular concentration 
could be attributed to the fact that the Laplacian is time reversal invariant (i.e. invari- 
ant under complex conjugation), so that the eigenf unctions are usually real- valued on the 
real submanifold M. Hence their complex zero sets are invariant under the (time reversal) 
involution a : (x, ^) —* {x, — ^) (the classical limit of complex conjugation). 

As a simple example of Theorem ll.il consider the circle S^. The geodesic flow is ergodic 
modulo the symmetry (x, ^) — >■ (x, — ^). The real eigenf unctions sin27r/cx, cos27rA;x are there- 
fore quantum ergodic. They complexify to the cylinder as sin 27Tkz, cos 2TTkz. The complex 
zero set of these holomorphic functions lies entirely on the set Q^^; = and become uniformly 
distributed with respect to 2d6 as A; — >■ oo. It will be checked in ^that the coefficient agrees 
with the result of Theorem 14.11 Note however that the complex eigenf unctions e^'^'^'^''^ are 
not quantum ergodic and have no complex zeros (see ^for further discussion). 

We now outline the main steps of the proof. As mentioned above, eigenfunctions ip\ 
of Laplcians of real analytic Riemannian manifolds admit holomorphic extensions to 
a maximal Grauert tube in the complexification Mc of M, which we will identify with a 
maximal ball bundle B*^M on which the adapted complex structure is defined. The square 

1^1^ of metric norm function \C,\g = \jYlTj=i 9^''{^)iiij is a smooth, strictly plurisubharmonic 
exhaustion function on B*^M. For < e < eo the sphere bundles S*M = dB*M are strictly 
pseudoconvex CR manifolds. We denote by 0{B*{M)) the class of holomorphic functions 
on this domain, and by 0{dB*[M)) the space of boundary values of holomorphic functions, 
i.e. the CR holomorphic functions. For each < e < eo, the restriction ^p'x\dB*{M) thus lies 
in the Hardy space 0^{dB*{M)) of square integrable CR functions. 
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A key object in the proof is the sequence of functions Ux{x,^) G C°°{B*M) defined by 

^7a(x,0:=#7T' {x,OeB:M, where 



(7) 



P\{X,0 ■= ll</'Al9%lJU2(9i?*^|^M) 



Thus, px{x,^) is the 'moving' L^-norm of ip'^ as it is restricted to the one-parameter family 
{dB*M} of strictly pseudo-convex CR manifolds. Ux is of course not holomorphic, but its 
restriction to each sphere bundle is CR holomorphic there, i.e. 

(8) ul = Ux\9B^MeO%dB:{M). 

Our first result gives an ergodicity property of holomorphic continuations of ergodic eigen- 
functions. 

Lemma 1.4. Assume that {(fix} 'is a quantum ergodic sequence of A-eigenfunctions on M in 
the sense of Then for each < e < eo, 

I^Ar-^ ^^^^^^^ leir^S weakly ^n L\B:M,u^). 

We note that uo^ = r^~^drdujdvol{x) in polar coordinates, so the right side indeed lies in 
L^. The actual limit function is otherwise irrelevant. The next step is to use a compactness 
argument from |SZj (see also |NVj ) to obtain strong convergence of the normalized logarithms 
of the sequence {|?7aP}- The first statement of the following lemma immediately implies the 
second. 

Lemma 1.5. Assume that \Ux\'^ ]i:{im)\^\g"''^^ ^ weakly in L\B*M,uj"^). Then: 

(1) ^ log |f/jp strongly m L^{B*M). 

(2) j-dd\og |f/jf ^ 0, weakly m V'{1, 1){B*M). 

Separating out the numerator and denominator of |f/j p, we obtain that 

(9) l5aiog|<^^P - ^ddlogpx^ ^ 0, (A,- ^ oo). 

The next lemma shows that the second term has a weak limit: 
Lemma 1.6. For < e < eo, 

j\ogpx{x,0^ L\b:M) as X^oo. 

Hence, 

^dd\ogpx^^dd\^\g^, (Xj^oo) weakly m V'{BIM). 

It follows that the left side of has the same limit, and that will complete the proof of 
Theorem 11.11 

The proofs of the lemmas are based on the properties of the analytic continuation of 
the wave kernel as a complex Fourier integral associated to the complexified exponential 
map |Bou| in52l IGLSj . Since (PxldB'^M = ^~^^E{i^)ipx, the analytic continuation of ipx 
is obtained by applying a complex Fourier integral operator of known order and symbol. 
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This allows us to connect growth and distribution of zeros in tubes to the dynamics of the 
geodesic flow. Although this paper treats only the ergodic case, it would be of interest to 
investigate complex zeros of A-eigenfunctions under other dynamical hypotheses such as 
complete integrability. It would also be interesting to investigate analogues for boundary 
value problems. 

Let us compare the results of this paper to earlier results of B. Shiffman-S. Zelditch |SZj 
and of S.Nonnemacher- A. Voros |NVj (see also [^) on complex zeros of eigenfunctions 
of ergodic quantum maps in Kahler phase spaces. These articles were concerned with the 
complex zeros of eigenfunctions of ergodic quantum maps acting on spaces H^{M^ L^) of 
sections of powers of a holomorphic line bundle over a a Kahler manifold (M, cu). The role 
of the Grauert tube was played by the disc bundle D* C L* of the dual line bundle, and the 
role of the CR manifolds S^M was played by the circle bundle X = dD*, which is a strictly 
pseudoconvex CR submanifold of L*. The norm function \^\g of {M,g) is thus analogous 
to the hermitian metric of (L, h), but the analogy is not very close. Indeed, the role of the 
geodesic role in the (M, g) setting was split in the Kahler setting between two dynamical 
systems: the action on the circle bundle X ^ M (which is the direct analogue of the 
geodesic flow in the Riemannian setting but is of course not ergodic); and an auxiliary ergodic 
symplectic transformation x of the Kahler manifold M. Ergodicity of joint eigenfunctions 
{snj} of the action and of the quantum map associated to x g^-^e ||sivj(2;)||^]v — * 1 in 
the weak sense (as in Theorem II. 4j) . where the norm is the pointwise norm relative to a 
hermitian metric on L with curvature equal to the Kahler form u. Weak convergence of the 
zero currents as in Lemma IT31 showed that j^Zg^. u, proving equidistribution of zeros 
relative to the Kahler form. In comparison. Theorem 11.11 shows that complex zeros of real 
ergodic eigenfunctions are not equidistributed relative to the analogous Kahler form Ug in 
the Riemannian setting, but rather to the relatively singular form stated in the theorem. 
The difference can be traced to Lemma 11.61 which indicates that the moving L^-norm px 
ends up playing the key role of the hermitian metric. 

1.2. A conjecture on ergodic real nodal hypersurfaces. We close the introduction 
by stating a conjecture on the real nodal hypersurfaces Z^^ = {x G M : ipj{x) = 0} on 
Riemannian manifolds with ergodic geodesic flow. We define the distribution of real zeros of 
ipx by integration 



with respect to (m — l)-dimensional Haussdorf measure dTi"^ ^ on the nodal hypersurface 
induced by the Riemannian metric of {M,g). 

Conjecture 1.7. Let {M,g) be a real analytic Riemannian manifold with ergodic geodesic 
flow, and let {fj} be the density one sequence of ergodic eigenfunctions. Then, 



More generally, we conjecture the same limit result for any quantum ergodic sequence of 
eigenfunctions. In the case of random spherical harmonics, the limit formula was proved 
in the PhD thesis of J. Neuheisel |Nej. We believe that (in a straightforward way) it can 



(10) 
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be extended to random orthonormal bases on any compact Riemannian manifold as defined 
above, which would give the asymptotic strengthening of |JLj in the real domain. 

There is of course a very wide gap between the equidistribution Conjecture 11.71 and the 
best known result on volumes (Q). Due to the singular concentration along the real zero set, 
it seems possible that Corollarv II .41 could have implications for the distribution of real zeros 
in the ergodic case. 
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in June 2005. We thank these institutions for their support. We would also like to thank 
Z. Rudnick, M. Sodin for suggesting improvements in the exposition and D. Jerison for 
encouragement to include Corollarv 11.31 

2. Analytic continuation to a Grauert tube 

In this section, we recall the relevant known results on analytic continuation of the wave 
kernel and A- eigenfunctions of a real analytic (M, g) to a Grauert tube. 

A real analytic manifold M always possesses a complexification Mc, i.e. a complex man- 
ifold of which M is a totally real submanifold (Bruhat- Whitney |BWj ). The germ of Mc 
along M is unique. In [0] , Grauert constructed plurisubharmonic exhaustion functions p on 
Mc, which define the Grauert tubes = {p < e} relative to p. 

In |GSll in52l ILSll ILS21 IGLSj . the complex geometry of Grauert tubes was brought into 
contact with the symplectic geometry of the cotangent bundle T*M and with the Riemannian 
geometry of real analytic metrics g. The key results (for the purposes of this article) are 
the following: First, a metric g determines a canonical plurisubharmonic function pg on Mq. 
It is defined on a maximal Grauert tube whose radius eo is often called the 'radius of the 
Grauert tube'. There is a symplectic diffeomorphism 

^p: M,^ B*M 

of the Grauert tubes with respect to pg to the ball bundles B*M C T*M with respect to g, 
which identifies pg with As pointed out in |LSlj (see also |GLSp . one may take ip''^ to 
be the complexified exponential map 

(x, e B*M ^ exp^. e M,. 

The map ip endows B*M with a complex structure Jg adapted to g. 

Proposition 2.1. |LS11 ITIC^ The adapted complex structure is uniquely characterized by 
the property that the complexified exponential map, 

{x, eB*M ^ exp^ e M, 

is a biholmorphism for < e < cq- 

The domains B*M are strictly pseudoconvex for e < eo, hence their boundaries S*M = 
dB*M are strictly pseudoconvex CR manifolds. We use the notation dB*M to emphasize 
the role of sphere bundles as boundaries of domains. The restrictions 

ifj, = exp(ie)-^ : dM, OB* 

of ip are then CR holomorphic diffeomorphisms. 
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As mentioned above, the metric norm function pulls back under ip^ to the function ^/p 
on Mc, which is known as the Monge- Ampere function |GSlj . It equals rc{z,z) where rc 
is the holomorphic extension of the distance function. In the cotangent picture, the metric 
norm function \^\g is smooth on B*^M\M and solves the homogeneous complex Monge- 
Ampere equation {dd\C,\g)"^ = there. In fact, the form dd\C,\g has rank m — 1 on B*^M\M, 
and its kernel is a smooth rank 1 sub-bundle of T[B*^M\M). The leaves of the associated 
('Monge- Ampere' or Riemann) foliation are the complex curves t + ir ^ '^'lit), where 7 is a 
geodesic, where r > and where T'y{t) denotes multiplication of the tangent vector to 7 by 
r. We refer to |LS1| for further discussion. 

2.1. Model examples. To better understand complexifications and Monge- Ampere func- 
tions, and in particular the limit form in Theorem ll.il we go over several (well-known) model 
examples. We note that these examples do not have ergodic geodesic flow, so they do not 
exemplify Theorem II. H but only the objects involved in it. 

(i) Complex tori: 

The complexification of the torus M = M™/Z"" is Mc = C""/Z"". The adapted complex 
structure to the flat metric on M is the standard (unique) complex structure on C". The 
complexified exponential map is exp^^{iC,) = z := x+i^, while the distance function r(x, y) = 
\x — y\ extends to rc{z,w) = \/ i^z — wY- Then ^(^z,z) = \/{z — zY = ±2i\'^z\ = ±2i\^\. 
Thus, the limit form is -dd\Qz\. 

(ii) E)"' |PW| IGSlj The unit sphere x1 + - ■ ■ + x^^j^^ = 1 in ]R"+^ is complexified as the complex 
quadric 

Si = {(zi, . . . , ^.) e C"+i ■.zl + --- + zl^, = 1}. 

If we write Zj = Xj + i^j, the equations become \x\^ — |^|^ = 1, {x, = 0- The geodesic flow 
G*(x, ^) = (costx + sint^, — sintx + cost^) induces the exponential map exp^C, = (cos|^|)x + 
(sin |.^|)^, which complexifies to 

expcx V -1^ = (cosh \^\)x + V-l(sinh 1^1)^- 
The distance function of S" of constant curvature 1 is given by: 

r(x, y) =2 sin"^ — - — = 2 sin"^(- a/ (x - y^), 
whose analytic continuation to x is the doubly-branched holomorphic function: 

rc{z, w) = 2 sin"^ ^ a/(z - wy. 

One branch gives the pluri-subharmonic function 

^{z) = rc{z,z) = 2sm-U\'^z\) = 2isinh-^ \Qz\) = zcosh"^ |z|^ {z G §^). 

Since 

exp£Vp(x,0 = cosh"^ |(cosh|^|)x + i(sinh|^|)|||2 
= cosh-^{(cosh - (sinh 1^1)2} 
= cosh-^cosh2|^| = 2|^|, 

the limit form is ^^SScosh^"*^ |zp on S^. 
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(iii) (See e.g. |KMj ). H" The hyperboloid model of hyperbolic space is the hypersurface 
in M."'^^ defined by 

Then, 

= {{z„ z^^,) G C"+^ : ^? + ■ ■ ■ 4 - 4+1 = -!}• 
In real coordinates Zj = Xj + , this is: 

{x,x)l-{^,0l = -1, {x,Ol = 

where (, )l is the Lorentz inner product of signature {n,l). The complexified exponential 
map is given by 

. II^IL . sin^ 

expcx(v -10 = cos(— + V-l( ii^ii^ K- 

Let 

We note that has two components according to the sign of ^Zn+i- The Monge- Ampere 
function is: 

^{z) = cos-\\\x\\l + m\l-rr)/V2. 

The radius of maximal Grauert tube is e = 1 or r = -n j Hence the limit form is 

^aacos-i(||x||i + ||e||i-vr)/v/2) on Ml. 

2.2. Analytic Continuation of the wave kernel. By the wave kernel of {M,g) we mean 
the kernel 

E{t,x,y) = J2('''^'^ji^)^jiy) 

of e**^. As discussed in |BSu1 OLB], the wave kernel at imaginary times admits a 
holomorphic extension to x M as 

oo 

(11) E(ze,C,2/) = 5^e-^^^Ci(C)^,(l/), {C,y)eM,xM. 

In the simplest (albeit non-compact) case of M", the wave kernel E{t, x, y) = J^„ e'^^^k'^^'^'-y^d^ 
analytically continues to t + ir, ( = x + ip & C+ x C" as the integral 

E{t + iT,x + ip, y)= [ e'(*+'")l«le'<«'"+'P-^>rfe, 



which converges absolutely for \p\ < r. At positive imaginary times and for (x, y) G M" x M", 
E{iT, X, y) is the Poisson kernel of the upper half space x M+, 



K{r, X, y) = ^1 + i^rj =t{t' + {x- y)')) " 
which visibly has a holomorphic continuation to ( = x + ip in the x variable for \p\ < t. 



11+1 

2 11+1 

2 
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On a general analytic Riemannian manifold, one has a similar integral formula for the 
wave kernel of the form 

(12) Eit,x,y)= [ e'*l^l''«e'<^'""P«"'(")U(t,a;,|/,0(^e 

Jt*m 

where \^\g^ is the metric norm function at x, and where A{t,x,y,^) is a polyhomogeneous 
amplitude of order 0. The analytic continuation of the wave group at imaginary times is the 
Poisson operator e~'^^, a Fourier integral operator with complex phase; for background on 
the microlocal analysis of such kernels, we refer to jTj (Chapter XI). The analytic continua- 
tion, 

(13) E{tTX,y) = / e--l«l«^e^<«'^^P«'(^)U(t,C,2/,0rfe iC = x + tp). 

of the Poisson kernel to the Grauert tube |C| < t defines a complex Fourier integral operator 
from L^(M) with values in holomorphic functions in where r < eg. Its canonical relation 
is the complexification of the canonical relation of the real wave group, i.e. 'graph' of the 
complexified geodesic flow at imaginary times 

(14) r,. = {(x, e, z, C) e T*M\0 X T*M,\0 : G'-{x, = {z, ().} 
We may (and will) also express the adjoint operator in the modified form 

(15) E*{iT,x,C)= [ e-"l«l«-e*<«'''"P-"'(^)U*(t,C,x,Of^^ {C = x + ip). 

Here, we use the phase — t|^|c,^ + i{^,eyip^^{Q) instead of —T\^\g^ — i{^,exp'^'^{x)). 

We can transport the complexified wave kernel to B*M x M using the complexified ex- 
ponential map (or ip): 

(16) ^(ze,(x,0,Z/) :=Eiie,exp^V^^,y). 

It is again a complex Fourier integral operator whose canonical relation Tir can be obtained 
from (fT^ by composing with ip. 

2.3. Szego projector. We denote by O'^'^^ {dB*M) the subspace of the Sobolev space 
W'^'^^ {dB*M) consisting of CR holomorphic functions, i.e. 

0'+'^{dB*M) = W'+'^{dB*M) n 0{dB*M). 

The inner product on 0^{dB*M) is with respect to the Liouville measure dfi^- There are 
similar spaces for dM^^ and composition with ip^ defines an isomorphism 

^* : O'^'^idB^M) 0"+V(aM,). 

We further denote by 

n, : L\dB*M) ^ 0\dB*M) 

the Szego projector for the tube B*M, i.e. the orthogonal projection onto boundary values 
of holomorphic functions in the tube. It is well-known (cf. |BoSj[ IMS| IGS2j ) that 11^ is 
a complex Fourier integral operator, whose real canonical relation is the graph As of the 
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identity map on the symplectic cone C T*{dB*M) spanned by the contact form a = ^-dx, 
i.e. 

= {(x,e;ra,), G dB*M), r > 0} C T*{dB*M). 

Alternatively it is a Toeplitz operator in the sense of Boutet de Monvel-Guillemin |BGj . The 
analogous Szego projector : L'^{dM^) (9°(9MJ is conjugate to 11^ in the sense that 

We now consider the restrictions o E{ie) from L'^{M) to 0{dB*M). Since is an 
]R+-bundle over dB*M, we can define the symplectic equivalence of cones: 

(17) : T*M ^ S„ i{x, = (a;, e^, 

The following result is the transport under -i/^^ to dB*M of results due to Boutet de Monvel 
(see also jOSS)- 

Theorem 2.2. [Boulin^ o E{ie) : L^{M) 0{dB*M) is a complex Fourier integral 
operator of order — associated to the canonical relation 

T = {{y,r],L,{y,r])}cT*MxE,. 

Moreover, for any s, 

n, o E{ie) : W'{M) ^ 0'+'^{dB*M) 
is a continuous isomorphism. 

2.4. Analytic continuation of eigenfunctions. We obtain the holomorphic extension of 
the eigenfunctions ipx by applying the complex Fourier integral operator Eiir): 

(18) E{it)^x = e-^Vv 

As usual, we can use the complexified exponential map to transport the ip^ to B*^{M): 

(19) ^f{x,0 = Vxiexp,V^O- 
By Theorem I2.2t we obtain: 

Corollary 2.3. |Bou| IGLSj Each eigenfunction (p\ has a holomorphic extension to B*M 
satisfying 

sup \<f^{x,0\<C,X"'+^e'\ 

{x,^)eB*M 

The fact that the holomorphic continuations of eigenfunctions can be obtained by applying 
a complex Fourier integral operator is the crucial link connecting the geodesic flow and the 
growth rate and zeros of <^a- 

2.5. Examples. We pause to consider some basic examples of holomorphic continuations of 
eigenfunctions. The simplest example is the flat torus ]R™'/Z™', where the real eigenfunctions 
are cos(fc, x), sm{k, x) with k G 2tt1I^. The complexified torus is 'U^ jTI^ and the complexified 
eigenfunctions are cos{k, (), sin(/c, () with ( = x + iC,. 

No such explicit examples exist in the ergodic case, but one can see the uniform analytic 
continuation of eigenfunctions very clearly in the case of compact hyberbolic quotients H^/F. 
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For simplicity we consider the two-dimensional case. Eigenf unctions can be then represented 
by Helgason's generalized Poisson integral formula (Hj . 

Jb 

Here, z G D (the unit disc), B = dD, and dT\ G V'{B) is the boundary value of ip\, taken 
in a weak sense along circles centered at the origin 0. Also, (2;, h) is the (signed) hyperbolic 
distance of the horocycle passing through z and 6 to 0. To analytically continue (p\ it suffices 
to analytically continue {z,h). Writing the latter as {C,b), we have: 

^^(C) = / e(^^+^)«-^)rfr,(6). 

Jb 

Using this representation, one could verify directly the estimates on growth of complexified 
eigenfunctions. 

3. The operators E{ie)*Il^aIl^E{ie) 

The following lemma will be used to reduce ergodicity properties and norm estimates of 
the complexified eigenfunctions to properties of the initial real eigenfunctions (px- We 
recall that the E, U notation refers to the B*M setting while E, U refers to the setting. 

Lemma 3.1. Let a e S^{T*M - 0). Then for allQ <t<tQ, we have: 

E{ie)*Ii,aIi,E{ie) E ^""^(M), 

with principal symbol equal to a(x,^) ^ . 

Proof. We observe that E{ie)*IleaIleE{ie) is a complex Fourier integral operator on L?'{M) 
associated to the canonical relation, 

F* o As. o r = {(y, r/, X, 0} C T*M x T*M : .,(x, = ie{y, v)}- 

Since is a conic symplectic isomorphism, it follows that F* o oF = At*mxt*m, i-e. that 
E{ie)*IleaIleE{ie) is a pseuododifferential operator. It follows from Theorem 12.21 that 

E{ie)*Il,aIl,E{ie) : W'{M) W'+'^iM), 

is a continuous linear map and hence that the order is — 

The principal symbol of E{ie)*IleaIleE{ie) equals a(x, ^) times the principal symbol of 
E{ie)*Il^E{ie). Indeed, by the calculus and Egorov theorem for complex Fourier integral 
operators, it equals the the principal symbol of E{ie)*Ilt,E{ie) times the translate of a under 
the canonical relation underlying Il^E[it). As discussed above, this relation is the symplectic 
identification of T*M - = S,. 

Thus, it suffices to show that the principal symbol of E{ie)*IleE{ie) equals \C,\g ^ . Since 
it equals the principal symbol of E{it)*Il^E{it) transported under ip, we can (and will) do 
the computation in the setting. The calculation could be done using the calculus of 
complex Fourier integral operators, but that would require a digression on symbols of Szego 
projectors and on the composition of the symbols of the three factors. It seems quicker to 
calculate the symbol from scratch by applying stationary phase for complex phase functions 
with only real critical points. 
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The calculation is well illustrated by the simplest case of the Euclidean wave kernel on 
M™. We then have: 

This is a complex Fourier integral operator with phase 

(20) ^oix, y, e, C, r) = -t{\^i\ + |6l) + ^(6, Xi + tp - y) - 2(^1, x - {x - ip)). 
The dxi integral produces 5(^1 + ^2), and then the d^i integral gives 

~ 2r'^ ^„ 1^1-"^ sinh2r|e|e-2-l«le-^<«'"-^>de 
In the last step we used the asymptotics of the inner integral 

(21) 

~ 2r^|er'^sinh2r|e|e-2-l«l 

Thus, E{iT)*IlrE{iT) is a Fourier multiplier by a polyhomogeneous function with leading 
term |4| 2 . 

We now show that the same result holds on a general compact analytic Riemannian man- 
ifold by reducing to the Euclidean case. Using the analytic continuation of the parametrix 
(fT!?|l . we have 



E{iTyUrE{iT){x,y) = /^..^^ /^.j^,^ Xj*^^ e-^(l«il«--+l«2l9,,)e-i(5i,cxp, '(C)>gi(6,exp/{C)> 

(22) 

X A* {t, C, X, ^i)A{t, C, y, i2)diidi2dV{z)duj 

The phase 



(23) vl> = _^(|^^|^-+ |^2|,,,) -«(ei,exp-i(C)) +*(6,exp,i(C)) 

is complex but has only real critical points given by the non-degenerate critical manifold 
- T*M: 

(24) Cq, = {{x,y,^i,^2,C = z + iTuj) : y = X, = -^2, z = 0, u = --^}. 

lul 

The Lagrange immersion 

: Cq, ^ T*M X T*M, ^^(x, y, ^1, ^2, C = -2 + ^™) = {x, 4*, y, -dy"^) 

is then real valued and, as mentioned above, we may apply stationary phase for complex 
phase functions with only real critical points (cf. |Hoj . Vol. 1) to evaluate the principal 
symbol. 

We recall that the principal symbol is the transport to ^t'MxT'M by cq, of the 1/2- 
density Cq a/ ^c* where (Tq is the principal term of the amplitude restricted to and where 
dc^ is the Leray density on C$ induced by the map iq, and the coordinate volume density 
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d$,id^2dV{z)du!. Above, ctq is the principal term of A*(t, (, x, C,i)A(t, (, y, ^2) 
Leray density is given by 



d, 



1 on The 





-1 




-1 


dxdC^i = 


dxd^i 


D(x,6,S/,C2,C) 





where we regard {x,C,i) as coordinates on We write ( = xi + ip in Riemannian normal 
coordinates based at x. Expanding first along the Dy rows and then along the "^'^^ columns, 
we obtain 



(25) 



D(p) 



dxd^i. 



'^'^^ column, it suffices to show that ^1^'^^^^ 



and that ^^'^^ 



set. In the calculation of '^'^^^i ^^e critical set we may put ^1 



Indeed, simple calculations show that the Dy derivatives only act non-trivially on , and 
that only has non-trivial derivatives under D^^ , D^^ . To eliminate the Dy rows and the 

Id = "if'i c on the critical 

and then the calculation reduces to Dl_^{{^, exp~^(xi — ip) — exp~^(a;i + ip)))\^^^Q = 0, which 
holds since exp~^ is the identity map when ( = x + ip is expressed in Riemannian normal 
coordinates at x. In calculating '^'^^^ = D'^^^{\C,2\gy +i{C.2,exp~^{()))\^=y^^^=_^^, the first term 
has a zero Hessian since all metric coefficients vanish in normal coordinates. The second 
is easily seen to be a constant multiple of the identity operator (the multiple is the same 
as in the Euclidean case). Finally, ^'^^^3 = "^-Daji exp~^(xi + ip)|^j=o,y=x = CId. Taking 
the determinant of the result gives Finally, we calculate this last determinant after 

restricting to the critical set, obtaining a constant multiple of 



det Dlp{^,exp/{ip))\j^^ 



CdetDl^{^,p)\ 



_ € . 



The last determinant is precisely the one that arises as the Hessian determinant in the 
stationary phase formula in (jUJ). The derivatives are taken on 5*™"^, hence we obtain 1^1™""^ 
times the normalized determinant, which is invariant under rotations and hence constant. 
Raising to the power — ^ completes the calculation. 

□ 



4. Proof of Lemma [HH and Lemma fTTT)l 

We now use Lemma 13.11 to reduce the quantum ergodicity and norm properties of the 
complexified eigenfunctions to properties of the original real eigenfunctions. 

4.1. Proof of Lemma 11.41 We begin by proving a weak limit formula for the CR holo- 
morphic functions u\ defined in (jH)) for fixed e. For notational simplicity, we drop the tilde 
notation although we work in the B*M setting. 

Lemma 4.1. Assume that {fx} is a quantum ergodic sequence. Then for each < e < cq, 

1 

l2,{dB*M)' 
That IS, for any a G C{dB*M), 

1 



weakly in L^{dB*M,dfXe 



dB*M 



a{x,^)dn^. 
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Proof. It suffices to consider a G C'^{dB*M). We then consider the Toephtz operator nealle 
on 0%dB*M). We have, 

(26) 

= e^^^^' ll<^AllL2%BjAf)(^(^e)*n,an,E(ze)(^j,(^j-)L2(M). 



m—l 



By Lemma 13.11 E{ie)*Il^aIl^E{ie) is a pseudodifferential operator on M of order 

m — l 

with principal symbol a\C,\g ^ , where a is the (degree 0) homogeneous extension of a to 
T*M — 0. The normalizing factor e^^'^^ I Iv^a I lL2(a_B*A/) same form with a = 1. Hence, 

the expression on the right side of may be written as 

^27) {E{ieyU,aU,E{ie)^j,(fj)L^M) 

{E{ie)*U,E{ie)Lpj,(fj)L2(^M) 

By the standard quantum ergodicity result on compact Riemannian manifolds with ergodic 
geodesic flow (see [ Shn| IZ2| [Z| ICVj for proofs and references) we have 

{E{ieyU,aU,E{ie)ipj,ipj)L2^M) ^ 1 f 
' ' (E(ze)*n,E(ze)¥P„¥P,)L2(M) y.,[dBlM) 

JdB*M 



m — l 



More precisely, the numerator is asymptotic to the right side times A 2 , while the denomi- 
nator has the same asymptotics when a is replaced by 1. We also use that ^ (ob'M) IdB*M '^^/^e 
equals the analogous average of a over dBi (see the discussion around (jHOI))- Taking the ratio 
produces (j2Hl)- 

Combining ^ and the fact that 

{U^aU^Upu'j) = / a\uj\'^dfi^ 

JdB*M 

completes the proof of the lemma. 

□ 

We now complete the proof of Lemma fl.4| i.e. we prove that 
(29) [ a\Ux\'uj^ ^ I 



B'M f^i{S*M) Jb*M 

for any a G C{B*M). It is only necessary to relate the Liouville measures dfir © to the 
symplectic volume measure. One may write dfir = ^\t=rXt^'^, where Xt is the characteristic 
function of B^M = {\^\g < t}. By homogeneity of Hr{dB;M) = r'"" Vi(<95*M). If 
a G C{B*), then J^,j^^auj"^ = jQ{jQ^,j^,jadfir}dr. By Lemma l^!Tl we have 

(30) 

/b|M «I^a|'^'" = loUoB^^M a\ul\^dfi,}dr ^ loijlAkf) loBfM adfirjdr 



Ibh 



ar ' LU 



Mi(9-B*Af) JB*M 

TT, |2 _ ,^ , 



^*-lim,_|t/,p = -^|e|-+i 
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4.2. Proof of Lemma II. 6L We actually prove the stronger result that 

\^ogpx{x,^) ^ \^\g^, uniformly in B^M as A ^ oo. 
A 

We state the weaker form because that is what we need in the proof of Theorem 11.11 
Again we drop the tilde notation for simplicity. 

Proof. Again using E{it)ip\ = e~^'^(f'^, we have: 

(31) = e^^^n,Eite)ipx,U,E{te)ipx)LHOB:M) 
Hence, 

(32) ^ log PA (a;, = 2|el,. + ^ log(E(ze)*n,E(ze)<^A, <^a). 
To complete the proof, we observe that 

(33) j\og{E{ieyil^E{ie)(px,(px) < C ^^^^-^ uniformly in e, 

where C is a constant independent of (e. A). Indeed, by Lemma l3.ll E{ie)*IleE{ie) is a 
pseudodifferential operator of order — ^^^^ with principal symbol |^| ~ . To obtain a uniform 
bound on {E{ie)*IltE{ie)ipx, ^Px) in any of the standard bounds for the norm of a pseudo- 
differential operator in terms of derivatives of the complete symbol would suffice. 

To take one such bound with a convenient reference, we write {E{ie)*IleE{ie)ipx, fx) as 
{l + \^)'^{{I-A)-'^E{ieyU,E{ie)ipx,ipx). Put A, := {I - A)-"^ E{ie)*U,E{ie). Since 
G vI/^(™~^), we may apply the Schur- Young bound of [H] (Vol. Ill, Theorem 18.1.11) to 
obtain 

(34) {Aefx, fx) < \\ AW L^^L^ < (sup |a,(x,OMM, 

where is the complete symbol of A^ relative to some choice of quantization a(x, D) of 
symbols. The complete symbol of A^ may be obtained by applying ((/ — A)~^^ to the 
representation in ()22|). It is clear that the complete symbol is smooth in the parameter e, 
hence the right side of ()34|) has uniform bound in e, proving ()34|) and therefore the lemma. 

□ 



5. Proof of Lemma [131 and Theorem 11.11 
The remaining step in the proof of Theorem 11.11 is the proof of Lemma 11.51 
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5.1. Proof of Lemma II. 5L 

Proof. The proof is similar to that of Lemma 1.4 of |SZj . We wish to prove that 

^- ■= -L log |f/^f ^ in L\B*M). 
Xj 

We argue by contradiction. If the conclusion is not true, then there exists a subsquence ipji^ 
satisfying W^jJIl^b^m) > S > 0. 

To obtain a contradiction, we first observe that ipj is quasi-plurisubharmonic (QPSH) on 
B*M, i.e. may be locally written as the sum of a plurisubharmonic function vj and a smooth 
function pj] equivalently iddipj is locally bounded below by a negative smooth (1, 1) form. 
Indeed we put 

:= ^log|(/?5^p, pj := -Pa,. 

We use the following fact about subharmonic functions (see |Ho[ Theorem 4.1.9]): 

Let {vj} be a sequence of subharmonic functions in an open set X C which have a 
uniform upper bound on any compact set. Then either Vj — ^ — oo uniformly on every compact 
set, or else there exists a subsequence Vj^ which is convergent in Lj^^{X). 

Since the proof is local, it also holds for open sets in manifolds, and in particular for 
X = B*M. 

We now verify that the hypotheses are satisfied in our example: 

• (i) the functions vj are uniformly bounded above on B*M] 

• (ii) limsupj-^^t;^ < 2|^|g. 

It suffices to prove these statements on each surface dB*M with uniform constants inde- 
pendent of e. On the surface dB*M, Uj = u^j. By the Sobolev inequality in {dB*M), 
we have 



Taking the logarithm, dividing by Xj, and combining with the limit formula of Lemma fl.6l 
proves (i) - (ii). 

We now settle the dichotomy above by proving that the sequence {vj} does not tend 
uniformly to — oo on compact sets. That would imply that ipj —oo uniformly on the 
spheres OB* for each e < cq. Hence, for each e, there would exist K > Q such that for 
k>K, 

(35) J_log|«^.^(z)| <-l. 

However, implies that 

\uj^{z)\ < e-"^^^" yzedB*M, 

which is inconsistent with the hypothesis that |m^-^(z)| — > 1 in T>'{dB*M). 

Therefore, the second half of the dichotomy holds, i.e. there must exist a subsequence, 
which we continue to denote by {vj,,}, which converges in L^{B*^^) to some v G L^(i?*^). By 
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passing if necessary to a further subsequence, we may assume that {vj,^} converges pointwise 
almost everywhere to v in B*^. Then, 

v{z) = hmsupwj^ < 2|^|g (a.e) . 

Now let 

v*{z) := limsupt'(w) < 

UI— »2 

be the upper-semicontinuous regularization of v. Then v* is plurisubharmonic on B*M and 
V* = V almost everywhere. 

Put ip* '■= V* — 2|^|g. Then ip* < 0, and the assumption ||'?/'jfe||Li(BjM) > 5 > implies that 

Us := {C e B^M : ^{0 <-5/2} 
has positive volume. Since ipj^ — >■ ip* in L^{Us), one has by jHoj Theorem 4.1.9 (b) that 

(36) limsupV^jJc/^ < V'ic/a < -S/2. 

fc— >oo 

Hence, there exists a positive integer K such that V'ife(C) ^ —S/2 for ( G Us, k > K; i.e., 

(37) |^,,(C)I <e-''^S C^Us, k>K. 

This again contradicts the weak convergence to 1. 

Therefore ||?/'j,.||Li(B*iv/) > S > leads to a contradiction, and the Lemma is proved. 



□ 



To complete the proof of Theorem 11.11 it suffices to combine the results that 
' -dd\og\Ux\' = T [Zx] - 7;^dd\ogM\\l,(SM.) ^ (weakly). 



and that (by Lemma fl.(i|l the second term tends to ^dd\^\g. □ 

5.2. Final remarks, (i) We check the numerical details in the case of the circle. 

The zeros of sm2'n'kz in the cylinder C/Z all lie on the real axis at the points z = 
Thus, there are 2k real zeros, and the Poincare-Lelong formula gives 

limfc^oo 2^98 log I sin 27rA;p = hmfc_,oo i Zlnti 

On the other hand. 



matching the other expression. 

As mentioned in the introduction, the complex eigenfunctions e^'^*'^^ have no complex 
zeros, hence Theorem II. II is false for them. The reason is that they are not quantum ergodic 
but rather localize on just one of the two components of the unit tangent bundle (the one 
with the same sign as k). Running through the previous calculation shows that the limit 
zero current for these eigenfunctions is ^dd$^ = rather than 
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(ii) One can obtain other formulae for the distribution of zeros in the ergodic case using 
the fact that the maps t + y/^s — > exp^^j) 57 (t) are holomorphic curves 7c (t + v^— Is) 
relative to the adapted complex structure for each geodesic 7. If one pulls back complexified 
eigenfunctions under 7c, then one obtains a holomorphic function in a strip around the 
real-axis. Its complex zeros are discrete and correspond to the intersection points Z^c fl 7^. 
Its real zeros are the intersection points of the real geodesic 7 with the nodal hypersurface. 
In connection with Conjecture ll.7t it is natural to conjecture that these intersection points 
become uniformly distributed on (M, g) when 7 is a uniformly distributed geodesic. 

(iii) We can give a simpler form to ^dd\^\g in dimensions m > 2. Let p{x,$,) = We 
note (with (HSl) that 

ddf{p) = f'{p)ddp + r{p)dpAdp, 

and that 

dp A dp = idp A a, ddp = —iug. 
It follows that in dimensions m > 2, we have 

fQQ\ ^ r- ^ dp A a 

(''^ ^^^v^=2V^"^+4V^- 
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